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We apply a recently developed quasiparticle self-consistent GW method (QSGW) to Gd, GdAs, 
GdN and ErAs. We show that QSGW combines advantages separately found in conventional GW 
and LDA+U theory, in a simple and fully ab initio way. QSGW reproduces the experimental 
occupied 4/ levels well, though unoccupied levels are systematically overestimated. Properties of the 
Fermi surface responsible for electronic properties are in good agreement with available experimental 
data. GdN is predicted to be very near a critical point of a first-order metal-insulator transition. 

PACS numbers: 71.15-m,71.10-w,71.20-Eh 



Most 4/ compounds belong to a class of materials 
whose electronic structure can be approximately de- 
scribed in terms of the coexistence of two subsystems — 
a localized / subsystem, and an itinerant spd subsystem. 
States near the Fermi energy E-p predominantly con- 
sist of the latter; 4/ electrons largely play a passive 
role except to spin-polarize the spd subsystem through 
an indirect exchange mechanism. Describing both sub- 
systems in the framework of ab initio electronic struc- 
ture methods, however, poses a rather formidable chal- 
lenge. The most widely used method, the local density 
approximation (LDA), has been an immensely success- 
ful tool that reasonably predicts ground-state properties 
of weakly correlated systems. The LDA is much less 
successful at predicting optical properties of such sys- 
tems, and its failures become serious when correlations 
become moderately strong. It fails catastrophically for 
open /-shell systems, leaving / electrons at Ep. To 
surmount these failures, a variety of strategies to ex- 
tend the LDA have been developed. These include exact 
exchange (EXX) [l[ , self-interaction-correction Q (SIC), 
LDA+U [1, i, [f|, and more recently LDA+DMFT (dy- 
namical mean-field theory) As a consequence, they 
have serious problems, both formal and practical. SIC, 
LDA+U and LDA+DMFT add nonlocal potentials to 
certain localized electrons in a special manner, leaving 
some ambiguity about how a localized electron state is 
defined, and how the double-counting term should be 
subtracted. Further, such approaches are specialized; 
they cannot remedy the LDA's inadequate description 
of itinerant spd subsystems (e.g. its well-known under- 
estimate of semiconductor bandgaps), which is the rele- 
vant one for transport properties. Thus, they are prob- 
lematic for 4/ materials such as the rare earth monop- 
nictides we study here. The standard GW (i.e. 1-shot 
GW as perturbation to the LDA, or G LDA W LDA ) sig- 
nificantly improves on the LDA's description of itinerant 
spd subsystems, but it has many shortcomings Q; and it 
fails qualitatively in open / systems, in much the same 
way as the LDA fails [7j. 

In short, the present status consists of an unsatisfac- 
tory patchwork of methods, each with successes in im- 



proving some property in one or another class of materi- 
als. Here we show, for the first time, that a recently devel- 
oped quasiparticle self-consistent GW method 
|ll| (QSGW) can reliably describe open /-shell systems, 
including the spd subsystem. The GW approximation 
is a prescription for mapping the non-interacting Green 
function to the dressed one, G° — > G. Formally, G can 
be calculated from any G°. QSGW gives a prescription 
to specify a (nearly) optimal mapping G — > G°, so that 
G° — > G — » G° — > ... can be iterated to self-consistency. 
At self-consistency the quasiparticle energies of G° co- 
incide with those of G. Thus QSGW is a self-consistent 
perturbation theory, where the self-consistency condition 
is constructed to minimize the size of the perturbation. 
QSGW is parameter- free, independent of basis set and of 
the LDA . It contains LDA+C7 kinds of effects, but no 
subsystem is singled out for specialized treatment; there 
are no ambiguities in double-counting terms, or in what 
is included and what is left out of the theory. We showed 
that QSGW reliably describes a wide range of spd sys- 
tems [9( . Its success in describing / systems is important 
because it is not known whether the GW method can 
reasonably describe correlated / electrons at all. 

In addition, this work brings out some new fundamen- 
tal points. First, the position of unoccupied / levels is 
systematically overestimated. Second, we predict that 
GdN is on the cusp of a new kind of first-order metal- 
insulator transition (MIT). Last, we show that the shifts 
in / levels that QSGW determines relative to LDA, lie 
outside the degrees of freedom inherent in the standard 
LDA+U method. The method enables us to reconstruct 
parameters used in the LDA+U theory, as we will show. 

Ref. [TlJ gives some formal justification as to why 
QSGW should be preferred to conventional self- 
consistent GW. The orbital basis and our development 
of the all-electron GW are described in Ref. Local 
orbitals (e.g. 5/ states) are essential for reliable descrip- 
tion of these systems in QSGW. It is also important not 
to assume time-reversal symmetry in the open / systems 

We considered the following 4/ systems: Gd, Er, EuN, 
GdN, ErAs, YbN, and GdAs. Gd and Er are metals, 
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FIG. 1: The QSGW energy band structure of Gd, GdAs, 
GdN, and ErAs. Right panels show DOS, together with 
experimental XPS and BIS data@, [U, [lj, [li| (circles). 
i?F=0eV. Colors indicate the spin character of the band (blue 
for majority and red for minority). Lattice constants were 
taken to be 5.00A(GdN), 5.80A(GdAs), 5.73A(ErAs) in the 
NaCl structure, and 3.64A(Gd) in hep. 



while the rest are narrow-gap insulators or semimetals. 
QSGW always shifts 4/ levels away from Ep. The elec- 
tronic structure around Ep is dominated by spd elec- 
trons, which we will consider later. Fig.Q] shows en- 
ergy bands and QP density of states (DOS), for some 
cases, together with XPS (X-ray photoemission) and BIS 
(bremsstrahlung isochromat) spectroscopies. 

/ subsystem: In all cases, stable ferromagnetic solu- 
tions were found with the 4/ element in the 3+ state: 
that is 6, 7, 11, and 13 / levels are occupied in Eu, 
Gd, Er, and Yb, respectively; the remainder are unoc- 
cupied. (Antiferromagnetic solutions were also found, 
but FM solutions are presented here to compare with 
Shubnikov-de Haas (SdH) experiments.) Occupied 4/ 
levels were always dispersionless, as Fig.Q] shows, while 
unoccupied states show some dispersion, reflecting their 
hybridization with the spd subsystem. Gd is the only 
4/ element for which the LDA and G LDA W LDA do not 
put / states at E?gfl8j|. because the majority 4/t states 
are filled and minority 4/^ states are empty. The two 
are separated by an exchange splitting (U-J in LDA+Z7 
terminology). For the remaining 4/ elements either the 
4/ T or the Af l level is partially filled. QSGW predicts 
large exchange splittings within this channel (controlled 
by different combinations of U and J); see e.g. the ErAs 
DOS in Fig.[T] Occupied / levels are generally in reason- 



TABLE I: QSGW 7 spin and orbital moments (bohr), average 
position of 4/ levels relative to Ep (eV), and corresponding 
peaks in XPS and BIS data (where available). When /"*• or 
f" states are split between occupied and unoccupied levels, 
average positions for both occupied and unoccupied are given 
(top and bottom numbers). When the occupied or unoccupied 
part of or f> levels consist of multiple states split about the 
average (see, e.g. ErAs bands in Fig.Q]), the range of splitting 
is denoted in parentheses. 
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able agreement with available XPS data (see TableflJ. In 
the two Er compounds (Er and ErAs), occupied f* and 
levels are fairly well separated. The shallower 4/^ 
levels likely correspond to the XPS peak between —4.5 
and — 5eV, and the 4/t levels to the broad XPS peak 
between —8 and — 10 eV shown in Fig.Q] In YbN, the 
separation between occupied /* and is small, and the 
XPS peak (whose width is ~3eV) probably corresponds 
to some average of them. More precise identification is 
not possible because multiplct effects are not included. 
In contrast, the unoccupied 4/ levels are systematically 
higher than observed BIS peaks, typically by ~3-4 eV. 
The only exception is ErAs, where the overestimate is 
closer to 1 eV. (This may well be an artifact of final-state 
effects in ErAs, as suggested in Ref. 0.) 

Overall, the 4/ T -4/ i splitting is 16.2 eV in Gd, and 
~18 eV in GdN and GdAs. This change in the splitting 
is reflected in the BIS-XPS data (12.2 eV for Gd, 14 eV 
for GdN and GdAs). The carrier concentration at Ep 
is larger in Gd than in GdAs, which results in a larger 
dielectric response, and more strongly screened U. 

The orbital moments folllow what is expected from 
Hund's rule. The spin moments are a little overestimated 
in the metals Er and Gd, following the trend observed in 
3c? magnetic systems such as MnAs [t| . 

Deviations from experiment can be qualitatively ex- 
plained as follows. QSGW 7 overestimates unoccupied 
states in sp semiconductors by ^0.2eV[§]. The over- 
estimate is somewhat larger in itinerant TM oxides such 
as SrTi0 3 and Ti0 2 (<leVV and is larger still (>leV) 
in the correlated oxide NiO [8j, [ll | . This can be under- 
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stood as a neglect of electron-hole interactions (excitonic 
effects). Short-range, atomic- like excitations shift peaks 
in Ime(o;) to lower frequency, increase the screening and 
reduce the strength of W . It is to be expected that the 
more localized states are, the stronger electron-hole cor- 
relations will be (4/ > localized 3d > itinerant 3d > sp) . 
Indeed, when Ime(w) is calculated through the Bethe- 
Salpeter equation, it is in dramatically better agreement 
with experiment in sp systems (see, e.g. Ref. |l9jh The 
expected corrections to W are consistent with the ob- 
served trends, both in the overestimate of unoccupied 
QPEs (increasing with localization), and the overesti- 
mate of magnetic moments. However, is also likely that 
vertex corrections to GW play an increasingly greater 
role as localization increases. It is beyond the scope of 
our ability to include either kind of vertex correction in 
QSGW at present. 

The spd subsystem comprises the states at Ep, which 
control electronic transport properties. TableHTl presents 
two of the de Haas-van Alphen (dHvA) frequencies ob- 
served in Gd. By comparing them to the calculated ones 
as a function of Ep, we can determine the shift in Ep 
required to match the dHvA data [2(| , and thus assess 
the error in those bands at Ep. TablclITl shows that 

the QSGW 7I and al should be shifted by 0.2 eV 

and — O.leV, respectively, consistent with precision of 
QSGW for itinerant systems @. The QSGW DOS at 
Ep (1.84 states/eV-atom) is slightly overestimated (1.57 
states/eV-atom plj). 

ErAs and GdAs may be viewed as slightly negative- 
gap insulators (semimetals), with an electron pocket at 
X compensated by a hole pocket at T (see Fig.QJ. Several 
experiments address the spd subsystem near Ep: 

(i) The As-p-like Tis band dispersion between T and 
X in ErAs. As seen in Fig.[TJ the X point has states at 
-1.2 eV and -0.8 eV, and split-off bands at -3.52 eV. A 
single band from T to X of width ~1.5eV was observed 
by photoemission (Fig. 7 of Ref. |j). 

(ii) SdH frequencies / and cyclotron masses m* (see 
Tablelllip. SdH measurements in Er n ssScn 39 As [22 1 
agree reasonably well with QSGW calculations for ErAs. 
Allen et al. [23[ estimated m*=0.17 from dc field measure- 
ments for ErScAs, which is consistent with the QSGW 
values m*(e j 4 T ) =0.16, m*(e A < )=0.13. Nakanishi et 
al. [iil identified two branches in the [100] direction from 
dHvA measurements in GdAs: m*=0.2 (/=246 T) for 



TABLE III: QSGW cyclotron masses m* , in units of the 
free electron mass m, and frequencies / (Tesla) for GdAs 
and ErAs. Three bands cross Ep near T (see Fig.[T}: the 
heavy hole (hi), light hole (h2) and split-off hole (sh). El- 
lipsoids at X have two inequivalent axes, e|j C and e^. Also 
shown are SdH frequencies measured for Ern b«Sco 39 As [22I]. 
Petukhov et al. showed that Sc doping has a modest effect 
on the Fermi surface, at least within the LDA(core-like 4/) 
approximation [fj . shT and shj. are not distinguished in ex- 
periments. Notation follows Ref. [53] except that T and J. are 
exchanged. 
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m*(e A ), and m*=0.26 (/=439 T) for m*(h 2 ). These 
are in good agreement with QSGW values in Table liTil 
Koyama et al. [25[ obtained m*=0.48 from a broad peak 
in cyclotron experiments. This may be understood as 
some kind of average of masses in Table lHll 

(Hi) The electron concentration in the pocket at X 
is controlled by the (negative) gap between T and X. 
The QSGW result agrees with experiment to within 
the reliability of the calculated bandgap (<~Q.2 eV): 
ErAs: 3.5 xlO 20 cm" 3 QSGW; 3 ± 1 expt, Ref [24] 
GdAs: 3.3 xlO 20 cm~ 3 QSGW; 2.3 expt, Ref [26] 
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TABLE II: Gd dHvA Frequencies (T) for a magnetic field 
oriented along the [0001] direction. 7I originates from the 
majority 6s band and al from the majority 5d band. Both 
lie in the FKM plane; they are depicted in Ref. [2p| ]. 
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FIG. 2: Energy bands near E F (E F =0) in GdN and ScN. 
For GdN, majority (minority) bands are drawn in dark 
(light) lines, (a) LDA+U for GdN ((7=6.7 eV, J=0.69eV) 
and LDA for ScN. A semimetal is predicted in both cases. 
(b) QSGW bands: both materials have positive gaps, (c) 
'scaled E' bands (see text). The ScN T-r and F-X gaps fall 
close to measured values (E C (X)-E V (T)=0.9±0.1 eV; E C (X)- 
J5„(X)=2.15eV^|). Right panel: DOS in QSGW (dark 
lines) and 'scaled £' (light lines). Scaling reduces the 4/ - 
4/T splitting, but the spd subsystem is little affected by it. 
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GdN is qualitatively similar to GdAs and ErAs, but 
there is some confusion as to whether GdN is an in- 
sulator or semimetal. Wacher and Kaldis measured a 
large carrier concentration (1.9x 10 21 cm 3 ). But, direct 
measurements of resistivity indicate insulating behav- 
ior 28J. Remarkably QSGW predicts two kinds of sta- 



ble, self-consistent solutions near the observed lattice 
constant a=5.00: one is insulating with E g (a)=E c (X)- 
J5„(r)=+0.2eV, and dE g /da = 2.7eV/A. It is stable for I 4 
a>4.96A. The other is a semimetal with E g (a)—— 0.2 eV 
and dE g /da = 2.1eV/A. It is stable for a<4.99A and 
at larger lattice constants when a tetragonal or trigonal 
shear is applied. Thus, a range of structures is found for 
which two solutions coexist, and the MIT is first- order. 
It can be connected with discontinuous changes in the 
dielectric function: e(ui— >0) diverges in the semimetallic 
phase, and approaches a constant in the insulating phase. 

To correct for QSGW's tendency to overestimate gaps 
slightly, we adopt the 'scaled £' approach [IoT| : we take 
a linear combination of the LDA and the QSGW po- 
tentials, (1 - a)x QSGW+ ax LDA. We found that 
a=0.2 can accurately reproduce experimental band gaps 
for a wide variety of materials; see Ref. [H] for III-V 
and II-VI semiconductors. Energy bands with a=0.2 
are shown in Fig.[2jc). For comparison, we also show 
the bands of ScN, whose electronic structure is similar 
to GdN near E F . Scaling brings the ScN T-X and X-X 
gaps to within ~0.1 eV of experiment, consistent with our 
general experience. It is reasonable to expect that GdN 
will be similarly accurate. Panel (c) shows that for GdN, 
£ 9 «0.05eV. Thus GdN is right on the cusp of a MIT. 
The spin-averaged X-X gap (1.48 eV for majority, and 
0.46 eV for minority) is in close agreement with 0.98 eV 
measured in paramagnetic GdN [29( . 

Comparison between QSGW and LDA+U: In both 
methods, the one-particle effective potential is written 
as V cS = V LDA + AV. In QSGW, AV is rather gen- 
eral; in LDA+U, AV is added only for the onsite // 
block which is specified by U and J[30j. For comparison, 
we extract the onsite // block (QSGW [/-/]) from V eS 
given by QSGW, by Fourier transform techniques 11 1. 
The spd subsystem near E-p almost perfectly recovers 
the LDA-like bands of the LDA+J7 method, Fig.^a); 
the / subsystem remains essentially unchanged from the 
full QSGW calculation, Fig.QJc). Thus, the QSGW H 
may be thought of as the "ultimate" LDA+J7 Hamilto- 
nian, with Hubbard parameters connecting all orbitals 
between all sites in a basis-independent way. Signifi- 
cantly, no choice of the standard LDA+C/ parameters 
can reproduce the QSGW (/^,/ T ) levels simultaneously, 
because their average position is essentially determined 
from the LDA. 

This work was supported by ONR contract N00014-02- 
1-1025. We are also indebted to the Ira A. Fulton High 
Performance Computing Initiative. A. N. Chantis would 
like to thank A. G. Pctukhov for fruitful discussions. 
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